Incomplete Andreev reflection in a clean SFS junction 



Jerome Cayssol 1,2 and Gilles Montambaux 1 
(l)Laboratoire de Physique des Solides, Associe au CNRS, Universite Paris Sud, 91405 Orsay, France and 
(2)Laboratoire de Physique Theorique et Modeles Statistiques, 
Associe au CNRS, Universite Paris Sud, 91405 Orsay, France 

We study the Josephson effect in a ciean Superconductor-Ferromagnet-Superconductor junc- 
tion for arbitrarily large spin polarizations. The Andreev reflection at a clean Ferromagnet- 
Superconductor interface is incomplete, and Andreev channels with a large incidence angle are 
progressively suppressed with increasing exchange energy. As a result, the critical current exhibits 
oscillations as a function of the exchange energy and of the length of the ferromagnet and has a 
temperature dependence which deviates from the one predicted by the quasiclassical theory. 



Current understanding of the Superconductor- 
Ferromagnet-Superconductor (SFS) Josephson effect 
is limited to small spin polarizations. In the case of 
conventional superconductors, the Josephson current 
is due to the Andreev conversion of singlet Cooper 
pairs into correlated electrons and holes with opposite 
spins propagating coherently in the ferromagnetic metal. 
Applying the Eilenberger equations 2] to a clean multi- 
channel SFS junction, Buzdin et al. 3] have predicted 
that this non dissipative current oscillates as a function 
of both the exchange energy splitting E ex and the length 
d of the ferromagnet, because of the mismatch 2E ex /hvp 
between the spin-up and spin-down Fermi wavevectors. 
This quasiclassical result assumes that the Andreev 
reflection is complete, as it is fully justified for weakly 
spin-polarized ferromagnetic alloys E ex <C Ep, Ep being 
the Fermi energy. First experimental evidence for such 
oscillating critical current has recently been reported in 
Nb-Cu-Ni-Cu-Nb junctions 0. The so-called 7r-phase 
state of a SFS junction || has also been observed using 
diffusive weak ferromagnetic alloys such as Cui__Ni~ |y] 
or Pdi_.Ni, UMM 

In the new field of spintronics, devices with high spin 
polarization are used in order to manipulate spin po- 
larized currents. In the recently discovered half metals 
(HM), such as CrC-2 and Lao.7Sr .3Mn03, the current is 
completely spin polarized because one spin subband is in- 
sulating. Ferromagnetic elements Fe, Co, Ni, also exhibit 
quite large spin polarizations. Anticipating the interest 
for large spin polarizations, de Jong and Beenakker 
have shown that in this case the Andreev reflection is not 
complete at a clean Ferromagnet-Superconductor (FS) 
interface, in contrast to the case of a clean nonmagnetic 
Normal metal-Superconductor (NS) interface. Even in 
the absence of impurity scattering, normal reflection may 
occur because of the diagonal exchange potential barrier 
between the ferromagnet and the superconductor. This 
suppression of the Andreev reflection affects preferen- 
tially the channels with large transverse momentum. As 
a result, the sub-gap conductance of a ballistic FS contact 
decreases quasi-linearly as a function of the spin polar- 
ization rj — E ex /Ep from twice the normal state conduc- 



tance (rj = 0) to zero (rj = 1), because of the progressive 
suppression of the Andreev process. Using this principle, 
a point-contact Andreev reflection technique has been de- 
veloped in order to mesure directly the spin polarization 
of materials 0]>E3> such as Lao.7Sr .3Mn03, Cr02, 
NiFe, NiMnSb, which were not easily accessible by spin 
resolved tunneling spectroscopy [T^. A huge amount of 
theoretical efforts has been devoted to transport proper- 
ties in a nanoscale FS contact ^3 ; E-l ; E-l ; E3 while few 
studies have considered the thermodynamical properties 
of FS heterostructures [lf|,|_Sj- 

In this Letter, we address the physics of the Josephson 
effect in a clean multichannel SFS junction in the range 
of arbitrarily large spin polarization. We show how the 
Josephson current is modified by the ordinary reflection 
induced by the ferromagnet in the crossover from a SNS 
(n = 0) to a S/HM/S junction (rj > 1). With increas- 
ing exchange energy, the Andreev reflection is suppressed 
for electrons propagating with a large incidence, so that 
the number of channels contributing to the total current 
decreases. This reduction of the number of "Andreev 
active channels" has furthermore a subtle effect on the 
Josephson current: although the FS conductance is al- 
ways reduced when ij increases [lf1 |. the critical current 
has a non-monotonic behavior, depending on the current- 
phase relationship of the suppressed channels. For large 
spin polarizations, the oscillations of the critical current 
depend separately on the product kpd and on the spin po- 
larization n. They are reduced and shifted with respect to 
the predictions of the quasiclassical theory Q in which 
only a single parameter, 2E ex d/ '(Tivf) — r/kpd, is rele- 
vant. For small spin polarizations, we naturally recover 
the quasiclassical results. In the HM limit E ex — ► Ep, 
the critical current vanishes because the Andreev reflec- 
tion is totally suppressed for all the transverse channels. 
In addition, we study the temperature dependence of the 
critical current for different values of the spin polariza- 
tion and of the length d of the ferromagnet. 

We consider a clean short SFS junction with a large 
number M of transverse channels and with a length d of 
the ferromagnetic region much smaller than the coher- 
ence length of the superconductor £ Q = hvp/A , where 



2 



A Q is the T = superconducting gap. The itinerant fer- 
romagnetism is described within the Stoner model by an 
effective one body potential V a (x) — —aE ex which de- 
pends on the spin direction, characterized by a = ±1. 
In the superconducting leads, V a (x) = 0. The supercon- 
ducting pair potential is A(x) =| A | e %x l 2 in the left lead 
and A(x) =| A | e~ %x l 2 in the right lead. In the absence 
of spin-flip scattering, the Bogoliubov-de Gennes equa- 
tions split in two sets of independent equations for the 
spin channels (u-\,v±) and (u|,i>f) 



H + V a (x) A(x) 
A(x)* -H* a + V*{x) 



= e(x) 



(1) 

where e(%) is the quasiparticle energy mesured from the 
Fermi energy |2jj . The kinetic part of the Hamiltonian 
H a = [(—ihd/dx — qA(x)) 2 — Ep]/2m, with the effective 
mass of electron and hole m, is expressed in terms of the 
vector potential A(x), which is responsible for the phase 
difference \ between the leads, and Ep — h 2 k F /2m is 
the Fermi energy. The Fermi velocities are identical in 
both superconductors and in the paramagnetic metal. 

Because both the pair and the disorder potential are 
identically zero in the ferromagnet, the eigenvectors of 
Eq. are electrons and holes with plane wave spatial 
dependencies. For a given transverse channel, the elec- 
tron and hole longitudinal wavevectors, k nr7 and h n — a 
respectively, satisfy 
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E n — Ep — e — er E ex , 



(2) 



where E n is the transverse energy of the channel. One 
may label the transverse channels by an angle 9 n which 
is the incidence angle of the corresponding quasiparticle 
trajectory 
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E n — 
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F sin 2 I 



> n = E F sin 2 
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From Eq. J5J), one sees that an electron with incidence 
9 n cannot form an Andreev bound state with a hole if 
E n = Ep sin 2 9 n > Ep ~ E ex . Therefore the electron is 
normally reflected as an electron with the same spin for 
angle 9 n > 9 n = arccos y/rj. Such a process is insensitive 
to the superconducting phase and thus carries no Joseph- 
son current. In the opposite case 6 n « 6„, the Andreev 
reflection is complete and supports a finite current. In 
the following, the former kind of channel is referred to as 
" Andreev inactive" and the latter as " Andreev active" . 

Recently, we have performed detailed studies of the 
spectrum of a single channel SFS junction for arbitrarily 
large exchange energies [2l|. Solving the Bogoliubov-de 
Gennes equations, the spectrum is found to be strongly 



modified in comparison to the quasiclassical spectrum 
[22l| because gaps open at \ — and \ — 7r - However, 
due to a cancellation between the corrections associated 
to each anticrossing, the current is almost unaffected up 
to very large spin polarizations 77 rs 0.95. The region in 
which Andreev reflection and ordinary reflection coexist 
is extremely small. As a result, the Josephson current 
through a single channel SFS junction is given to great 
accuracy by the formula for perfect Andreev reflection 
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i(x,k F d,r],9 n = 0) = — ^ 
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for r\ < 1 and it is zero for 77 > 1. The parameter a = 
(yl + 77 — sfl — r/)kpd is the phase shift accumulated 
between an electron and a hole located at the Fermi level 
during their propagation on a length d. 

In the present paper, we generalize this result to trans- 
verse channels with finite angle 9 n , in the more realistic 
case of a finite width SFS junction. The crossover be- 
tween Andreev active and inactive channels occurs in 
a narrow window of incidences in the vicinity of 9 V = 
arccos y/rj. Below this cut-off, the current carried by a 
single Andreev active channel is 
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and it is zero for 9^ > arccos y/rj. In order to treat large 
exchange splitting, one has to take into account the exact 
band structure (here a simple isotropic parabolic band) 
and to express the phase shift between an electron and 
its Andreev reflected hole by 



kpd cos 9 r , 
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instead of using the linearized form 



rjkpd 2E ex d 



cos n hvpeost 



. (6) 
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The transverse channels considered above are indepen- 
dent because V a (x) is translationaly invariant in the 
transverse directions. Thus, the total current is the sum 
of the currents carried by each of them. As we assume 
a large number of channels, the discrete sum over n can 
be replaced by an integral over the angle 0. Calculating 
the total current, one has to restrict the integration over 
Andreev active levels only, so that the angular integral 
has to be limited by the upper cut-off n = arccos y/rj 



I(X, k F d, 77) = 



kpS 
~2~^~ 



dO sin cos i(x, kpd, 77, 



(8) 



FIG. 1: Current-phase relationships at zero temperature for 
a — 7r/4 obtained for several pairs (77, kpd). In the quasiclas- 
sical approximation, the current is a function of the single 
parameter a and does not decrease with increasing rj. The 
current is given in units of I = 7rAo/(ei?jv) 

where S is the cross section area of the ferromagnet. 

This expression, together with Eqs. (J5J|BJ|, is the cen- 
tral result of this Letter. It gives the Josephson current 
I(x, kpd, 77) of a clean SFS junction in the regime of arbi- 
trarily large spin polarization. Examples of current-phase 
relationships are shown in Fig. ^ In the limit of small 
polarization rj = E ex /Ep — ► 0, we recover the quasi- 
classical current-phase relationship [3| in which all the 
transverse channels contribute because 6 V — > 7r/2. 
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FIG. 2: Zero temperature critical current I c {f]) as a function 
ofn = E ex /Ep for different lengths of the ferromagnet, kpd = 
1, 5, 10. The current is given in units of I — 7rAo/(e_Rjv). 

Increasing the spin polarization, we study how the crit- 
ical current evolves from the case of a weakly spin polar- 
ized junction to the S/HM/S junction. As shown in Fig. 
13 the critical current has a non trivial oscillatory behav- 
ior as a function of exchange splitting for a given length, 
namely for fixed kpd. The number of oscillations oc- 
curing during the crossover from the SNS (77 = 0) to the 
S/HM/S junction (77 = 1) decreases when kpd is lowered. 
In the limit of an ultra-small junction kpd ps 1, there are 
no oscillations because the phase shift in Eq. © tends to 
zero, and all transverse channels carry the same SNS cur- 
rent with maximal value i Q = 27rA/0o, where <p = h/e 
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FIG. 3: Zero temperature critical current I c (n) as a function 
of kpd (thick lines), for different values of the spin polariza- 
tion 77. As 77 increases, the exact current deviates from the 
quasiclassical estimate (dashed lines). The current is given in 
units of I = 7rA/(ei?jv). 

is the flux quantum. Consequently, the reduction of the 
total current is only governed by the upper cut-off in Eq. 

m 

h = Mi (i - v ) = -^-(1 - ??). (9) 

This linear reduction of the current with increasing the 
exchange field is quite reminiscent of the almost linear re- 
duction obtained in Ref. for the conductance of a FS 
nanocontact |2^| . The total number of transport channels 
M = kpS/Air is large and determines the small normal 
state resistance Rn — h/(2e 2 M) of the heterojunction. 
The natural unit for the critical current is I = ttA/cRn, 
namely the one of a short clean SNS junction. 

Fig. |3| represents the critical current as a function of 
the length d of the ferromagnetic region, for different spin 
polarizations. We find that the oscillations are reduced 
and shifted with respect to the quasiclassical calcula- 
tion. There are two reasons for these deviations. Firstly, 
trajectories with large incidence are progressively sup- 
pressed. Secondly, the phase shift between electrons and 
holes for a given channel [Eq. JHJl] depends on the partic- 
ular band structure and differs from the linearized version 
a n = r\kpdj cos9 n . For large d, the oscillations decay 
slowly at zero temperature. In real situations, they are 
expected to be severely reduced when d exceeds the ther- 
mal length hp = Tivf/T or the phase coherence length 

We finally consider the effect of a finite temperature 
on the critical current. We have adopted the BCS tem- 
perature dependence of the order parameter A(T) = 
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FIG. 4: a) Critical current as a function of the spin polar- 
ization r\ at T = 0.9 T c . It vanishes for particular values of 
the spin polarization, when the junction undergoes a — ir 
transition. Io(T) = 7rA(T) 2 /(4e_RjvT c ) is the critical cur- 
rent for a SNS junction, b) Critical current [in units of 
Iq — 7rA /(ei?jv)/ as a function of the reduced temperature 
T/T c for values of n corresponding to the maxima of Fig (a), 
c) Critical current as a function ofT/T c for different values of 
n corresponding to the — n transitions. All curves correspond 
to a short junction with fc^d = 10. 

A tanh(1.74 v /T c /T - 1), and the exchange energy is as- 
sumed to be temperature independent. For T » T c , Fig. 
2Ja) shows that the critical current oscillates with the 
spin polarization 77 and cancels out for some values of 
77. In this temperature range, the current-phase relation- 
ship is sinusoidal I(x) = Ic sin% and the current vanishes 
identically when I c is zero. These cancellations are asso- 
ciated to transitions between the 0-phase state and the 
7r-state of the junction. For fixed parameters fcp-d and 
77, the critical current decreases monotonously with in- 
creasing temperature T, as shown in Figs. 0fb,c). This 
temperature dependence is very sensitive to the spin po- 
larization. For polarizations corresponding to — tt tran- 
sitions, I C (T) decreases exponentially with temperature 
[Fig. Efc)] whereas a much more slower decrease is ob- 
tained for the local maxima of the critical current [Fig. 
12b)]. 

We have studied the Josephson current of a clean SFS 
junction for arbitrary large spin polarizations. The two 
physical effects involved are the reduction of the num- 
ber of active levels participating in the Andreev process 
and the use of the non-linearized band structure. In any 
experiment with strong ferromagnetic elements or nearly 
half metallic compounds, the critical current oscillations 
should be affected by these effects. Firstly, the oscilla- 



tions depend separately on the spin polarization 77 and 
on the product kpd instead of the single combination 
rjkpd as suggested by the quasiclassical theory. Secondly, 
when the temperature is increased from zero to the criti- 
cal temperature, the local minima of the current are more 
strongly suppressed than the local maxima. 

We thank Igor Zutic for useful comments. 
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